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Electron and proton heating by solar wind turbulence 
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Abstract. Previous formulations of heating and transport associated with strong mag- 
netohydrodynamic (MHD) turbulence are generalized to incorporate separate internal 
energy equations for electrons and protons. Electron heat conduction is included. En- 
ergy is supplied by turbulent heating that affects both electrons and protons, and is ex- 
changed between them via collisions. Comparison to available Ulysses data shows that 
a reasonable accounting for the data is provided when (i) the energy exchange timescale 
is very long and (ii) the deposition of heat due to turbulence is divided, with 60% go- 
ing to proton heating and 40% into electron heating. Heat conduction, determined here 
by an empirical fit, plays a major role in describing the electron data. 
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1. Introduction 

The solar wind displays a highly non-adiabatic temper- 
ature profile, requiring some process(es) to provide addi- 
tional heat sources. One possible, and successful, source of 
heating comes from magnetohydrodynamic (MHD) turbu- 
lence present in the solar wind [Coleman, 1968]. An ac- 
tive MHD turbulent cascade [MacBride et al, 2008; Marino 
et al. , 2008] transfers energy from the large-scale fluctuations 
down to small scales where kinetic processes dissipate the en- 
ergy as heat. Previous theories [Zhou and Matthaeus , 1990; 
Marsch and Tu, 1993; Oughton and Matthaeus, 1995; Zank 
et al, 1996; Matthaeus et al, 1999; Smith et al, 2001, 2006; 
Matthaeus et al, 2004; Breech et al, 2005; Isenberg et al, 
2003; Breech et al, 2008], have been able to account for ra- 
dial evolution of fluctuation level, correlation scale, cross he- 
licity and temperature in a specified background solar wind 
flow. These theories (with some exceptions [e.g., Cranmer 
et al, 2007]) have focused only on proton temperature, ig- 
noring heat conduction and, in effect, assuming that dissipa- 
tion of turbulence occurs only through proton kinetic chan- 
nels. While neglect of proton heat conduction is justified, it 
is not obvious why all turbulent heating should impact pro- 
tons alone, nor is it clear why, or whether, proton-electron 
energy exchange can be neglected. 

Electrons provide additional heating by carrying the bulk 
of the solar wind heat flux and through collisions with the 
protons. Turbulence models often neglect electrons in favor 
of the protons as the protons help set the scales of interest. 
This is particularly clear in the case of the momentum con- 
tent and the mass density of the solar wind plasma. However 
the internal energy content of the electrons is not negligible 
compared to that of the protons. 

For this work, we seek to understand how electrons and 
protons conspire to heat the solar wind. Two possible av- 
enues immediately open up to explore these issues: 1) an 
empirically based approach to compute the heating rates 
and 2) a modelling approach based on turbulence theory. 
A companion paper [Cranmer et al, 2009] follows the first 
avenue, e.g., examines effects of electrons through a more 
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empirically based approach. Here we focus on modelling 
the heating of the solar wind through MHD turbulence the- 
ory. For the first time, effects of electron heat conduction 
are included in a turbulence transport model that extends 
beyond the inner heliosphere. This provides an important 
step toward realism and completeness in the turbulent heat- 
ing models, and enables the use of additional observational 
constraints. 

Under the assumptions of a spherically expanding, radi- 
ally constant solar wind, the equations for the steady-state 
temperatures of electrons, T e , and protons, T p , may be writ- 
ten as 
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where r is heliocentric distance, U is the bulk solar wind 
speed, ni e:P \ is the electron (proton) number density, Q{ e , P } 
represents turbulent heating per unit mass, q e is the elec- 
tron heat flux vector (proton heat flux vector has been ne- 
glected), and ks is Boltzmann's constant. The T e —T p terms 
in Eqs. (l)-(2) model Coulomb collisions taking place over 
a timescale r [e.g., Priest, 1982]. 

In previous models [Zank et al, 1996; Matthaeus et al, 
1999; Smith et al, 2001; Breech et al, 2008] the equation for 
evolution of the temperature has been supplemented by an 
equation for the turbulence energy, an equation for the cor- 
relation or energy-containing scale and an equation for the 
cross helicity. These will be revisited further below. Here we 
will focus on issues surrounding the inclusion of the separate 
temperature equations for protons and electrons. In partic- 
ular, three questions arise by writing down these equations: 
(1) How does the electron heat flux vary with distance? (2) 
How much turbulent dissipation goes into heating the elec- 
trons versus the protons? and (3) Over what timescale do 
the electrons and protons experience energy exchange cou- 
plings, such as Coulomb collisions, that tend to equilibrate 
their temperatures? 

In this paper, we explore the physical issues surrounding 
these questions within the context of a turbulence transport 
model for the solar wind. The physical issues may be further 
complicated by the effects of pickup protons and latitudinal 
variations present in the solar wind. As a result, we limit 
ourselves to the high latitude fast wind, which displays less 
dependence on latitude, at least above 35 degrees or so [Mc- 
Comas et al, 2000]. We also limit ourselves to examining 
these issues primarily in the inner parts of the heliosphere 
(r < 10 AU), which avoids complications due to pickup pro- 
tons. In future work, we will examine how pickup protons 
and latitudinal variations affect the key issues listed above. 
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2. Electron Effects 
2.1. Electron Heat Flux 

One of the new effects we include is the electron heat flux. 
Generally speaking the heat flux associated with the proton 
distribution is regarded as unimportant [Braginskii, 1965]. 
For the case of the solar wind, a simple way to see this is to 
assume that the proton heat flux is due to a beam with frac- 
tional number density, /beam, moving at the Alfven speed, 
Va, relative to the bulk proton population. This leads to a 
heat flux on the order of /beamVj. If we take the scale of 
its divergence to be R, the heliocentric distance, then the 
relative contribution of the divergence of the proton heat 
flux to the proton temperature given in Eq. (2) is of the 
order /bcamUf/i?. We compare this result to the heating 
term (see §2.3 below), which is of order Z 3 /\, where Z is 
the fluctuation amplitude and A the correlation scale. Since 
Z ~ Va/2 and A ~ 7?/100, while /beam cannot be greater 
than unity and is probably < 1/10, it is clear that for pro- 
tons the heat flux is much less important than the heating. 
Electrons, on the other hand, with their high thermal speed 
and anisotropic distribution functions, are affected greatly 
by heat flux. Within the solar wind, the electron heat flux is 
primarily along the magnetic field. Heat flux transverse to 
the magnetic field is severely reduced as the magnetic field 
effectively acts as an insulating blanket. 

Modeling q„ — the electron heat flux along the magnetic 
field — requires approximating the solar wind as either col- 
lision dominated or collisionless. A collision dominated 
model, ^ = — «n Vm T e (r) [Spitzer and Harm, 1953] pro- 
duces temperatures that are too large at 1 AU compared to 
observations [Hollweg, 1976]. Collisionless models [Hollweg, 
1976] fare much better and may be adequate for most pur- 
poses. Nonetheless, the collisionless models can still suffer 
from missing non-local effects [Canullo et al, 1996; Scudder 
and Olbert, 1979]. 

Rather than employing a theoretically based model, here 
we adopt a more empirical approach in order to handle these 
difficulties in specifying . We derive from observa- 
tional data by applying a best fit procedure to the elec- 
tron heat fluxes given by Pilipp et al. [1990] for 0.3 AU to 
1AU and Scime et al. [1994] for 1AU to 5.4 AU. Figure 1 
shows the observations and the fitted curve. Introducing 
x = log 10 (r/i?i au), the final fit is given by 



log 10 q (r) = -2.3054 - 2.115a - 0.58604a; 2 



(3) 



and is assumed to be a reasonable approximation out to 
r ~ 100 AU. Note that past 10 AU, pickup proton effects 
may invalidate this assumption. 
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We remark that the use of this electron heat conduc- 
tion data may introduce some uncertainty into our results. 
In particular, the measurements shown in Figure 1 were 
adopted from Scime et al. [1994]. The observations come 
from the initial Ulysses cruise phase from 1 to 5 AU, while 
Ulysses was in the ecliptic plane. To our knowledge, no 
velocity selection was performed on the observations. We 
therefore assume the measurements apply to all wind speed 
intervals and latitudes. Note that in some later observa- 
tions, Scime et al. [1999] found no significant variation in 
at higher latitudes, which are dominated by the fast wind. 
We will return to this point later when discussing our re- 
sults. 

2.2. Collision Effects 

Protons and electrons can exchange heat through 
Coulomb collisions, though they may do so over extremely 
long timescales. They each can also exchange energy with 
the electromagnetic field, and therefore indirectly with each 
other, through wave-particle interactions. These kinetic ef- 
fects may be thought of as producing effective collisions. 
The observational fact that protons and electrons frequently 
maintain different temperatures in the solar wind [Hund- 
hausen, 1972] provides clear evidence that this coupling is 
not very strong. Generally, the scales for observing colli- 
sional effects correspond to several AU, as has been found 
recently [Kasper et al., 2008] in a comparison of fast and 
slow wind characteristics where proton-ion collisions are the 
central consideration. 

Taking into account this background, we adopt a (some- 
what) arbitrary approach for setting the collision timescale, 
t. For this work, we set r as a constant equal to the plasma 
transit time to some distance (e.g., 10 AU). This choice pro- 
duces only a weak interaction between protons and electrons, 
but does allow for collisional energy exchange effects to be 
seen in the model results. 

2.3. Partitioning Turbulent Heating 

As the solar wind plasma evolves, MHD turbulence trans- 
fers energy from low wavenumbers to high wavenumbers, 
where the energy is dissipated. The dissipated energy heats 
both the protons and electrons. Since we are modeling the 
supply of energy from large scales to small scales, a pro- 
cess believed to be controlled mainly by large-scale MHD 
processes, we cannot, on the basis of this analysis, distin- 
guish the channel or sequence of kinetic processes that ab- 
sorbs the energy. Possibilities are kinetic Alfven waves [Lea- 
mon et al, 1998; Cranmer and van Ballegooijen, 2003; Bale 
et al, 2005; Gary and Borovsky, 2008], whistlers [Gary et al., 
2008], and nonlinear dissipation in current sheets [Sundkvist 
et al, 2007], to name a few. To account for the data, how- 
ever, we will have to employ a reasonably correct partition- 
ing of the cascaded energy into the heat functions Q p and 
Q e that appear in the temperature equations (1) and (2). 

Specifically, let Q be the total energy dissipated by the 
turbulence. We then define 



Qv — fp 



Qe = (l- f v )Q, 



(4) 



Figure 1. Observed electron heat flux (cgs units) [Pilipp 
et al, 1990; Scime et al, 1994] and fits to the data. 



where f p is the fraction that determines the amount of tur- 
bulent heating that goes into the protons (Q P ). The rest of 
the heat is given to the electrons (Q e ). We make no distinc- 
tion here between the core, halo, and strahl populations of 
the electrons as the turbulence heats all electrons. 

In principle, f p can be determined by the underlying ki- 
netic physics [Gary and Borovsky, 2004], which operates at 
scales far smaller than the energy-containing scale described 
by the turbulence model. f p may also change with increasing 
distance [see the campanion paper Cranmer et al., 2009, for 
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more details on how f p may change]. We neglect these fea- 
tures of the problem for now and simply adopt f p — const. 

Some evidence suggests [Leamon et ai, 1998, 1999] that 
f p ~ 0.60 at 1 AU, if one associates helicity-sensitive pro- 
cesses with proton cyclotron absorption, and the remaining 
heating as due to a process like Landau damping, which can 
affect both protons and electrons. There has also been the 
suggestion [Matthaeus et al, 2008] that a stronger cascade 
drives proton instabilities more robustly, leading to more 
absorption near proton cyclotron scales. However, lacking 
a firm quantitative theoretical basis, such empirically based 
arguments must be regarded with caution at present. There- 
fore we revert to a purely empirical choice of f p , to be sup- 
ported or contraindicated by the results of the model, shown 
below. 

For the total heating rate due to the cascade, Q, we em- 
ploy a von Kdrmdn and Howarth [1938] style phenomenol- 
ogy giving 

7 

Q ~ T (5) 

with Z 2 being (twice) the total energy available to the tur- 
bulence and A being the length scale of the largest turbulent 
structures. Both of these quantities must be extracted from 
the turbulence model. Note that this dissipation model is 
rather robust and has been successfully applied to very dif- 
ferent regimes, such as the super Alfvenic solar wind [e.g., 
Breech et al, 2008] and the corona [e.g., Dmitruk et al, 2001; 
Verdini and Velli, 2007]. 

3. Turbulence Model 

We adopt the following turbulence model [Zank et al, 
1996; Matthaeus et ai, 1999, 2004; Smith et al, 2001, 2006; 
Isenberg et at, 2003; Breech et al, 2005, 2008] which com- 
prises three equations describing the transport of turbulence 
quantities; one equation for (twice) the total turbulent en- 
ergy, Z 2 , 
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correlation (similarity) scale, A, 



d\_ +Z PE PI 
dr 11 U aUZ 2 ' 



(6) 



(7) 



of a = 0.8. For this work, we adopt lower values of a = 0.5 
based on the results of hydrodynamical simulations [Pear- 
son et al, 2004]. Even lower values, down to a = 1/8, may 
also be consistent with those simulations. 

The turbulence dissipates energy, which then goes into 
heating of the solar wind and manifests as elevated electron 
and proton temperatures. The total turbulence heating rate 
is given as 
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Figure 2. Solution of the turbulence model used in 
this paper. The boundary values are specified at 0.3 AU 
as Z 2 = 5000 (km/s) 2 , A = 0.03 AU, and a c = 0.6. 
Other parameters are set as U = 774 km/s, C s h= 0.25, 
a — 2/3 = 0.5, and Mao = —1/6. Pickup protons ef- 
fectively turn on near 10 AU. The solutions (solid lines) 
compare reasonably well to the observed Ulysses values 
from near solar minimum [Bavassano et ai, 2000a, b]. 
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Figure 3. Model solutions for the proton and electron 
temperatures computed without collisions and without 
(8) electron heat conduction. The proton solution matches 
the observations well, but the electron solution misses 
almost all of the observed data. 



/ K) = 



[(l + a c ) 1/2 ±(l-*c) 1/2 ], (9) 



and /'(c c ) — — / • Driving of the turbulence comes 

from shear, modeled through C s h, and pickup protons, Epi, 
with shear being more dominant in the inner heliosphere 
(say, r < 10 AU) and pickup protons more dominant in the 
outer heliosphere. M = 1/2 relates to the underlying tur- 
bulence geometries and a D — —1/3 approximates the nor- 
malized energy difference (kinetic minus magnetic) of the 
fluctuations. More details on the model itself, and its pa- 
rameters, can be found in Breech et al. [2008]. 

The rightmost term of Eq. (1), af + Z A /XU, controls the 
turbulent dissipation. The constant a (taken = 2/3 herein) 
modulates the efficiency of the dissipation, with higher val- 
ues resulting in a more efficient dissipation (e.g., more ac- 
tive turbulent cascade) . Previous use of the turbulent model 
[such as Breech et al, 2008; Smith et al, 2006] used values 
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Figure 4. Model solutions for the proton and electron 
temperatures computed with collisions but without elec- 
tron heat conduction. The proton solution is reason- 
able, but the electron solution falls below most of the 
observed data. Collisional effects are also apparent, with 
the two temperatures begining to equalize in the outer 
heliosphere, although the protons stay warmer. 
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which is partitioned between the electrons and protons ac- 
cording to Eq. (4). The factor of 1/2 appears because Z 2 is 
twice the energy density of the fluctuations. 

4. Model Results 

We can numerically solve the transport model, Eqs. (1)— 
(8), and use the results to compute the temperature values. 
Figure 2 shows the solutions of the turbulence model. The 
computed solutions compare well against Ulysses observa- 
tions of 1 hour cadence from solar minimum [Bavassano 
et al, 2000a, b]. The model can also produce other solu- 
tions, which bracket the observed values [see Breech et al., 
2008]. 

The transport solutions can then be used to compute the 
proton and electron temperatures using Eqs. (1) and (2) and 
taking f p as a constant 0.6. The temperature solutions given 
below use observations from the SWOOPS [Bame et al., 
1992] experiment on Ulysses corresponding to the same time 
as the observations shown in Fig. 2. The temperature ob- 
servations were made in the high latitude, fast wind during 
solar minimum conditions. For the proton temperatures, we 
have used the geometric mean of the two temperature values 
given in the data archives [see Cranmer et al, 2009, for more 
information]. For the electron temperatures, SWOOPS pro- 
vides temperature data for the core, halo, and total electron 
populations. We compare our results against the total elec- 
tron temperature as the turbulence heating should affect all 
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Figure 5. Model solutions for the proton and electron 
temperatures computed with both collisions and electron 
heat conduction turned on. The proton solution remains 
reasonable. The electron solution displays a "shelf-like" 
region between 1 and 10 AU. The solution matches the 
observations better than not including the heat conduc- 
tion (cf. Fig. 4). 
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Figure 6. The electron temperature solution from Fig- 
ure 5 showing the "shelf" region more clearly. The so- 
lution fits most of the data, although there are outliers, 
particularly around 3 AU. The "shelf" region is due to 
heat conduction (see text). 



electrons. For the model solutions, we take the initial values 
of T p = 2.0 x 10 6 K and T e = 4.0 x 10 5 K at 0.3 AU. 

Figure 3 displays the temperature solutions without heat 
conduction and almost without collisions. The collision 
timescale r was set equal to the time for plasma to transit 
to 100 AU, which effectively removes collision effects from 
the results. The proton solution agrees very well with the 
observed data, but the electron solution misses almost all 
the data. 

Figure 4 shows the temperature solutions without the 
electron heat conduction, but with collisions turning on at 
10 AU (i.e., r was set to the transit time to 10 AU). The 
proton solution mostly agrees with the observations, but 
does show some disagreement near 5AU where the proton 
temperature solution is slightly cooler than the observations 
indicate. Collisions evidently cause the slight disagreement 
as this disagreement is not present in Fig. 3. This may pro- 
vide evidence that setting the collision time to be the transit 
time to 10 AU is too low a value. 

In contrast to the proton temperature solution, the elec- 
tron solution becomes slightly better when collisions are in- 
cluded. The energy lost from the protons profits the elec- 
trons. Nonetheless, the electron solution still misses most of 
the observed data. 
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Figure 7. Model solutions for the proton and elec- 
tron temperatures computed with electron heat conduc- 
tion turned on, but without collisions. The shelf in the 
electron temperature solution is still present, as is the 
crossover of the electron and proton solutions. The initial 
value of the electron temperature was raised to 6.0 x 10 5 K 
to maintain the good agreement with observations. 
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Figure 8. Solutions for the electron temperature using 
different profiles of the heat flux vector, q\\ . The baseline 
solution uses the heat flux vector given in Eq. (3). The 
other two solutions use baseline heat flux multiplied by 
2 and divided by 2. Possible variations in the heat flux 
should be encompassed by those two extremes, the figure 
demonstrates that a higher heat flux will make the shelf 
region more pronounced, while a lower heat flux lessens 
the shelf. 
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Including both electron heat conduction, as determined 
from equation (3), and collisions produces interesting re- 
sults (Fig. 5). The electron temperature now features a 
"shelf" region between 1 and 10 AU where T e (r) does not 
decrease as rapidly as it does outside that region. Interest- 
ingly, the electron temperature solution crosses the proton 
temperature solution near 5AU, after which the electrons 
are actually warmer than the protons. Collisions then begin 
to take hold, resulting in higher proton temperatures and 
lower electron temperatures. 

Figure 6 shows only the electron solution from Figure 5. 
The shelf region matches the observed electron tempera- 
tures better than the solutions obtained when heat conduc- 
tion is neglected. While there is some variability in the ob- 
servations, the solution goes through most of the available 
data. Within this region, the heat conduction evidently al- 
lows some heat to "pile up" to allow the electrons to stay 
warmer than they may otherwise be. The shelf arises from 
the electron heat flux shifting from a form close to r~ 2 to a 
steeper form in the outer heliosphere. 

We note that the existence of the shelf region must be 
handled carefully. The shelf originates from the model solu- 
tions and compares favorably with the data we have. How- 
ever, no high latitude data is available beyond about 5 AU. 
Further data would be required to confirm the existence of 
shelf. Additionally, while we chose the electron tempera- 
ture data to match published data sets [Bavassano et al., 
2000a, b], it may be possible that there are data selection 
effects that lend to the appearance of the shelf. The tem- 
perature data itself was taken solar minimum conditions and 
has a cadence of 1 hour. Selecting data with other cadences 
or different solar cycle conditions could easily obscure the 
shelf [see, for instance, Phillips et al, 1995]. 

The shelf is due solely to electron heat conduction, as 
Fig. 7 reveals. For the solutions shown there, collisions 
were turned off (i.e., r is set to the transit time to 100 AU). 
The shelf is still present, although the initial value of the 
electron temperature had to be increased from 4.0 x 10 5 K 
(as in the other solutions) to 6.0 x 10 5 K to maintain similar 
agreement with the observations. Note that the proton so- 
lution displays better agreement with the observations due 
to the lack of collisions. 

The shelf is somewhat sensitive to the underlying electron 
heat flux vector, gy, used to compute the solutions. Fig- 
ure 8 shows the electron temperature solutions computed 
using three different profiles for qu . Equation (3) provides 
the baseline profile, which yields the electron temperature 
solution shown in Figure 6). The other two solutions use a 
heat flux vector where the baseline gy was multiplied by 2 
and divided by 2. Since Scime et al. [1999] found small varia- 
tions in the heat flux with latitude, any variation in the heat 
flux measurements should lie between these two extremes. 
The electron temperature solutions show that increasing the 
heat flux makes the shelf more pronounced. Lowering the 
heat flux lessens the shelf. 

5. Conclusions 

For this work, we have added the effects of electrons to 
a turbulence transport model for the solar wind. Elec- 
tron effects manifest in the division of turbulent heating 
between protons and electrons, collisions between protons 
and electrons, and heat conduction of the electrons. We 
find that adding an empirical model for electron heat con- 
duction makes a dramatic difference in the solutions. A 
"shelf-like" region appears between 1 and 10 AU where the 
electron temperatures do not decrease as rapidly as they do 
outside that region. The electrons stay warm enough to ac- 
tually become warmer than the protons. Observations of 
the electron temperatures seem to support the presence of 
the shelf. We also find that collisions between the protons 
and electrons may not be important until well past 10 AU. 
Allowing collisions to take effect closer in appears to cause 
discrepancies between the model's proton temperature solu- 
tion and observed data. 
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